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ABSTRACT 

We present a general star formation law where star formation rate depends upon effi- 
ciency a, timescale r of star formation, gas component a g of surface mass density and 
a real exponent n. A given exponent n determines r which however yields the corre- 
sponding star formation rate. Current nominal Schmidt exponent n s for our model is 
2 < n s < 3. Based on a gravitational instability parameter Qa and another dimen- 
sionless parameter fp = (P/Ga 2 ) 1 ^ 2 , where P = pressure, a c = column density of 
molecular clouds, we suggest a general equation for star formation rate which depends 
upon relative competence of the two parameters for various physical circumstances. 
We find that Qa emerges to be a better parameter for star formation scenario than 
Toomrc Q-parameter. Star formation rate in the solar neighbourhood is found to be 
in good agreement with values inferred from previous studies. Under closed box ap- 
proximation model, we obtain a relation between metallicity of gas and the efficiency 
of star formation. Our model calculations of metallicity in the solar neighbourhood 
agree with earlier estimates. We conclude that metallicity dispersion for stars of same 
age may result due to a change in efficiency through which different sample stars were 
processed. For no significant change of metallicity with age, we suggest that all sample 
stars were born with almost similar efficiency. 
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1 INTRODUCTION 

It was realized by Kenicutt (1989) that there occurs non- 
linear increase in the star formation rate near the threshold 
surface density corresponding to Q-parameter. This shows in 
fact why rate of star formation R is large in the spiral arms 
of several galaxies. For example, in M 51 and NGC 6946 
(Rydbeck, Hjalmarson & Rydbeck 1985; Lord 1987; Tacconi- 
Garman 1988) , gas densities in the spiral arms are larger by 
a factor of two indicating deviations in the usual power law 
exponent (n ^> 2) of Schmidt (1959, 1963). In fact, star for- 
mation in many spiral galaxies under extreme conditions of 
low gas density and low disk self-gravity present a challange 
to all current theories for disk star formation (Ferguson et 
al. 1996). New star formation laws have therefore been pro- 
posed (see e.g. Talbot & Arnett 1975; Dopita 1985; Wyse 
1986; Silk 1987; Koppen & Frohlich 1997). 

However, a global star formation law has been put to 
doubt (Hunter & Gallagher 1986). For a general star for- 
mation scenario one may refer to Zinnecker & Tscharnuter 

(1984) ; Zinnecker (1989 and references therein). Thus, many 
interesting variants on the simple star formation laws in- 
clude e.g. self-propagating star formation (stochastic) in 
Gerola & Seiden (1978); Seiden & Gerola (1982); Dopita 

(1985) ; Coziol (1996); self-propagating star formation in Ari- 



moto (1989) and Hensler & Burkert (1990a, 1990b); star 
formation bursts (stochastic) in Matteucci & Tossi (1985). 
Kriigel & Tutukov (1993) and Tutukov & Kriigel (1995) have 
used one-zone dynamical code without radial dependence of 
the variables to study the conditions for bursts of star forma- 
tion. In the later paper, using one-zone code, they studied 
other types of burst of star formation in a galactic nucleus 
different from periodic bursts. Further, surface gas density 
threshold for star formation has been discussed in Kenicutt 
(1989). 

Wyse & Silk (1989) have discussed an extended Schmidt 
model with R dependence on surface gas density a g and lo- 
cal angular frequency fi(r) for both atomic and molecular 
gases respectively with n = 1 and n — 2. Wang & Silk (1994) 
have presented recently a self-consistent model (considering 
total gas surface density) for global star formation based on 
the gravitational instability parameter Q < 1 by Toomre 
(1964). In the solar neighbourhood, the model agrees with 
(i) the observed star formation rate, (ii) the metallicity dis- 
tribution among G-dwarfs, and (iii) the age metallicity re- 
lation for F-dwarfs. The model may be compared for star 
formation rate in galactic disks with Schmidt law with an 
exponent of about 2. Star formation rate depends also on 
the epicyclic frequency. A natural cut-off for Q = 1 in the 
star formation rate results. However, their analysis is heav- 
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ily based on Q < 1 criterion which has been put to question 
in relation to non-radial instabilities in the galactic disks 
that may play more fundamental role when magnetic field 
supported by azimuthal gas motions (thus resulting thermal 
instability not at all related to Q) is taken into considera- 
tion (Elmegreen 1993) . We note (see e.g. Figures 4 and 6 of 
Wang & Silk 1994) that the star formation does proceed in 
the regions where Q > 1. A natural question to ask is: how 
does star formation occur when Q > 1 and consequently 
the system has attained the state of gravitational equilib- 
rium? We attempt here to precisely answer this question 
and provide a scenario to circumvent this natural cut-off in 
the star formation process (see Sect. 3 in the text for details). 
The fact that star formation occurs via gravitational insta- 
bility was also suggested by Fall & Efstathiou (1980). The 
Q-regulation near its threshold value has been discussed by 
Dopita (1985) and Silk (1992). Silk (1995) has argued that 
local self-regulation of star formation may help expain the 
initial mass function of stars and that global self-regulation 
can account for the rate of star formation. Effect of envi- 
ronment on the gas content and rotation curves of disk may 
play a crucial role in determining star formation rates and 
histories. 

Elmegreen (1995) has discussed critical column densi- 
ties for gravitational instabilities and for cooling to diffuse 
cloud temperatures. It has been shown that the fundamen- 
tal scale for star formation in the outer regions of galaxies 
(in spiral arms) and in the resonance rings are related to 
the local unstable length. Since critical gas density for grav- 
itational instability scales as local density, inner regions of 
galaxies have higher star formation rate beyond threshold 
density. 

The consideration of magnetic field changes velocity dis- 
persion by a factor of two pushing Q > 1, (i.e. stable region). 
Incorporating this with the fact that there occurs shear in- 
stability of magnetised gas in the azimuthal direction, one 
is stimulated to think that Q < 1 may not be the only crite- 
rion for cloud formation which leads to star formation. An 
alternative suggestion for cloud formation because of energy 
dissipation accompanied by shear instability thus leading to 
star formation (even if Q > 1) has been given (Elmegreen 
1991a, 1993, and below for details). Macroscopic thermal in- 
stabilities and various cloud formation mechanisms are re- 
viewed in Elmegreen (1991b). We assume that instability 
parameter suggested by Elmegreen (1993), i.e. Qa < 1 (in- 
stead of Q < 1) is the criterion which determines occurrence 
of significant cloud formation instabilities. A natural conse- 
quence of our analysis is that star formation proceeds in 
the regions where one has Q > 1. It may be noted that in 
these regions (system being gravitationally stable) an alto- 
gether different cloud formation mechanism (leading to star 
formation) as suggested by Elmegreen (1993) is asked for. 
We shall present subsequently the evidence in support of 
our assumption. The outline of the paper is as follows: we 
give a general law for star formation rate in Section 2. In 
Section 3, we suggest a general equation for star formation 
rate which depends upon two fundamental parameters Qa 
and fp (defined in the text). We also give a comparison of 
star formation rate in the solar neighbourhood and timescale 
of gas depletion. Variations in the star formation rate and 
metallicity distribution in the solar neighbourhood are dis- 



cussed in Section 4. The Section 5 gives discussions and a 
resume of our results. 



2 STAR FORMATION RATE 

We write the star formation rate in the form 



(1) 



where a — efficiency of star formation which also depends 
upon n, t = timescale of star formation, a g — surface den- 
sity of gas composed of atomic and molecular components, 
n — an exponent. Clearly, r _1 is related to growth rate of 
instability (Goldreich & Lynden- Bell 1965). A review of re- 
cent observations of the history of star formation and its rel- 
evance to galaxy formation and evolution has been discussed 
by Kennicutt (1996). For the evolution of the global star for- 
mation history measured from the Hubble Deep Field one 
may refer to Connolly et al. (1997). Gravitational instability 
of galactic disks has also been studied by Elmegreen (1979), 
Cowie (1981), Ikeuchi, Habe & Tanaka (1984) and Bizyaev 
(1997). However, gravitational instabilities in the presence 
of turbulence are discussed in Bonazzola et al. (1987) and 
Leorat, Passot & Pouquet (1990). It is found that super- 
sonic turbulence may be strong enough (in some cases) to 
hold the Jeans criterion for gravitational instability. As a 
result, it may stop gravitational collapse. In this scenario, 
star formation takes place in molecular cloud complexes at 
places where the turbulence evolves to subsonic phase. In 
the present analysis we do not aim to discuss the instability 
criteria and their relevance to star formation (which are cer- 
tainly interesting topics of research at present), instead we 
aim to obtain a general star formation law with small num- 
ber of adjustable parameters. We assume neither infall nor 
radial flow in the disk. We consider gravitational instability 
due to axisymmetric perturbations (for non-axisymmetric 
case, one may refer to Goldreich & Lynden-Bell 1965) with 
magnetic field in the azimuthal direction which gives rise to 
shear instability in a magnetised gas. Groth rate of instabil- 
ity is now expressed as 



Li. < — k V e ff — 2-KGogk + k 2 



(2) 



where k is the wave number and k is the epicyclic frequency, 
v e ff is the effective velocity dispersion for ambient Alfven 
speed such that 



Veff 



I 2 2 n1/2 

(V Jeff+VAlf) , 



(3) 



v being the velocity dispersion without magnetic field and 
fcj — u) c (l + s — 2r) 



7e// 



to + uj c {?> — s) 



(4) 



Here 7 = ratio of two specific heats, lo c = cooling rate (see 
e.g. Elmegreen 1993 for details). The parameter Q is written 
as Q = KV e f f /irGug . Gravitational instability requires that 
both Q < 1 and k are smaller than a critical value 



^f^i + O-eT' 2 ) 



(5) 



Due to thermal instability, if 7 e // reaches large negative 
values (such that v\f f < 0), it implies no critical (or min- 
imum) wavelength for gravitational perturbation in the ra- 
dial direction. This makes Q 2 < 0. However, we do have 
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a maximum wavelength of the perturbation. Thus, equa- 
tion (2) shows the absence of Q-threshold for azimuthal 
instability which means that all Q-values provide unstable 
growth. Q-threshold may appear only if j e ff(uj) becomes a 
constant. Therefore, for the present treatment, we demand 
Qa = 2\/2Av e ff/irGa g < 1 for growth of gravitational in- 
stability but we are well aware that thermal and shear in- 
stabilities (along azimuthal direction) are capable of deter- 
mining cloud formation leading to star formation even if 
Q > 1. This amounts to replacing k by 2\/2A in the original 
Q (a scale transformation for Keplerian disk). 
Maximum of J 2 occurs at 

k max = 2V2A/v ef fQ A (6) 

(A = Oort shear constant), which provides maximum growth 
rate as (Wang & Silk 1994) 

,1/2 

(7) 



2 V / 2A(1 



Qa) 



Qa 

Since, r ~ |aw x | _1 , one gets from equations (1) and (7) 
a(2V2A)"a^(l-Q A r /2 



R n = 



(8) 



Following Wang & Silk (1994), we define a function f c = 
a g /<j c , <Jc = column density of individual molecular clouds. 
But, however, the relation between individual cloud forma- 
tion and star formation is complicated. Even cloud forma- 
tion process is not well known. The assumption that star 
formation rate results due to gravitational instability natu- 
rally demands for its relation with cloud formation process. 
Elmegreen (1990) has shown that gravitational instabilities 
generally form giant molecular clouds faster than due to ran- 
dom collisions. Cloud formation followed by star formation 
in the interstellar medium is certainly not the purpose of 
our investigation. Under the assumption that only gravita- 
tional instability is predominant, small cloud collisions may 
lead to large molecular clouds wherein star formation ensues. 
It is then natural to think that within an order of magni- 
tude cloud formation timescale or equivalently cloud colli- 
sion timescale and growth rate of local instability timescale 
are similar. With this scenario, Wang & Silk (1994) derive 
the expression for collision time between two clouds. We, 
thus, make use of their result and write collision time be- 
tween two clouds as 



-l 

coll 



,{2y/2A) 



(9) 



i 

coll 



• u max , and one gets 



OcQa 

In view of the above, t 

Qa~(1-/ c 2 ) 1/2 . (10) 

It is to be noted that this may not reflect the general prop- 
erty of the interstallar medium, e.g. other types of instabili- 
ties namely thermal instability and Parker instability might 
also contribute and affect the timescale of star formation 
(subsequently other physical quantities). Making use of eq. 
(10) into eq. (8), star formation rate is now expressed as 



R 



a(2V2A) n aZf? 



(11) 



(1 -/»)»/* • 

Eventually, in this form eq. (11) now assumes the conver- 
sion from column density to density using the galactic scale 
height. Let us write equation (11) in the form 



Table 1. 


Variation of shear constant A with surface density 


Distance A 


log A 




logo's 




(km s hpc ^ ) 




(M QP c-' 2 ) 




1 


105 


2.0212 


100 


2.0000 


2 


30 


1.4771 


3 


0.4771 


3 


20.9 


1.3202 


5 


0.6990 


4 


19.7 


1.2945 


10 


1.0000 


5 


19.1 


1.2820 


10.5 


1.0212 


6 


18.2 


1.2601 


10.2 


1.0086 


7 


17.2 


1.2355 


10 


1.0000 


10 


13.8 


1.1399 


7 


0.8451 


12 


11.5 


1.0607 


5 


0.6990 


14 


9.6 


0.9823 


4 


0.6021 


16 


7.9 


0.8976 


3 


0.4771 


18 


6.5 


0.8129 


2 


0.3010 


20 


5.44 


0.7356 


1 


0.0000 



dlnR din A 
n— = n + n— h n 



d 



d In a „ 



d\no„ 



d In a„ 



hi 



fc 



(i - m i/2 



dlnR 



„ dlnA 

+ 



fc 



(12) 



d In o g d In a g 1 — /f 

where n s stands for the nominal Schmidt exponent. The sec- 
ond term in equation (12) appears because for spiral waves 
epicyclic frequency is expressed through 

K = Ko(o" 9 /o"()) 1/2 

and shear constant A is 



(13) 



A = A {2 - (Jg/ao). 



(14) 



Non-axisymmetric gravitational perturbation of a magnetic 
gaseous disk has been discussed by Elmegreen (1987) who 
has obtained eq. (13) and eq. (14). Here, A and a represent 
equilibrium values of the shear rate and the surface mass 
density (see also Waller & Hodge 1991). It is easy to see 
that for vanishing shear constant, equation (12) reduces to 
equation (19) of Wang & Silk (1994, hereafter WS). It may 
be regarded as generalised version of WS equation in the 
sense that there is an additional term on the right side which 
is certainly non-zero. We calculate the second term on the 
right of equation (12), i.e., d log A/d log a g ~ 0.54 using least 
squares method. Data reported in Table 1 have been taken 
from Einasto (1979) and WS. Since mostly f c is very small 
compared to unity (see e.g. Table 5) for present Galactic 
disk, we conclude that the nominal Schmidt exponent n s 
for our model corresponds to 2 < n s < 3 for the Galaxy. 
For usual Schmidt law n s lies between 1 and 2. The other 
normal spiral galaxies of the Galaxy type are supposd to 
follow the same signature. 



3 GENERAL EQUATION FOR STAR 
FORMATION 

We suggest, that two fundamental parameters (Elmegreen 
1993) which determine star formation may be put in the 
form 



R n = aa n + P(af P ) n 



(15) 
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where a = 2^/2 Ao 2 g /Q A o c and / P = (P/Go 2 ) 1/2 , dissipa- 
tion and P is the pressure. The dimensionless pressure fp 
(defined originally by Elmegreen 1993) is the square root of 
the ratio of cloud collision rate to the gravitational instabil- 
ity rate and so is a measure of the relative importance of 
cloud collisions. In our discussions presented here, we make 
use of some interesting results of Elmegreen (1993). Both, 
Qa and fp now determine star formation rate. When both 
are large (i.e., Qa > 1 and fp S> 1), either thermal instabil- 
ity (macroscopic) triggers star formation or cooling which 
is very effective reduces Qa till gravitational instabilities 
switch over. When both are small (Qa < 1 and fp <ti 1), 
gravitational instabilities form clouds quickly but star for- 
mation is hampered due to lack of energy dissipation. How- 
ever, when Qa > 1 and /p < 1, star formation proceeds via 
random cloud collisions triggered by thermal instability and 
rate R is determined by second term in equation (15). This 
is believed to occur at galactic radii r < 4 kpc and r > 8 kpc 
where one observes Q > 1 (for an observed Q-distribution in 
the Galaxy, see e.g. WS). When Qa < 1 and fp 3> 1, grav- 
itational instability is primarily responsible for both cloud 
and star formation at all radii. In this case, the first term in 
equation (15) determines R. It is found that at all radii star 
formation is governed by relative competence of either of 
these terms. It also becomes evident that virtually Qa < 1 
(or Q < 1) is not an absolute criterion for star formation, 
instead it proceeds continuously until requisite ingredients 
are fuelled in and physical conditions are met. In fact, one 
observes significant star formation even when Qa > 1 in 
the Galaxy. Thus, the process of star formation can be visu- 
alized through the equation (15). It may be remarked that 
our model does not take into account the galactic bulge com- 
ponent (Oort 1977) which might contribute to galactic gas 
dynamics in the inner region below 0.1 kpc. 
We write star formation rate as 



dt n ' 



(16) 



where 5 is the fraction of mass returned to interstellar 
medium from the stellar content. From equations (15) and 
(16) we get 

a(2V2A) n a 2n f3(2y/2A) n cr, 



V? — fp -{T=s) IF' (17) 



QaoS 

Assume that parameters fp and cr™ are independent of time. 
Let us write Equation (17) in the form 

[q(1 - S) n (2V2A) n + (3(1 - 8) n (2V2A) n fp] dt n 

(w 2 r /2 



fc 



df c . (18) 



Integrate equation (18) to obtain 



,n ,n /i _ r2\n/2 

— + —fp = sin-^/D + constant, (19) 

7~a Tfj Jc 



where we have put 

r- n = a(l-8) n {2V2A) n 
r7 n = f3(l-8)"(2V2A) n . 



(20) 



If we express f g = a g /ai, fd = Oi/oc, where the subscript 
i denotes initial values of quantities, we can write equation 
(19) as 



r2 r2 \n/2 

JgJd) — sin- 1 (f^f^) +CO nstant.(21) 



r r r = _0_ 

At t = 0, LT g — <7i (f g = 1), one obtains the value of constant 
in equation (21) as 



constant - 



(i - fid 



2 \n/2 



+ sin- 1 (/5). 



(22) 



Now, equation (21) becomes 



(i+/p)- 



(i- 



fn f n 
J 9 Jci 

fir' 2 



+ sin 



J ci 

_1 (/3) 



-sin-^/S) 



(23) 



where we have put r™ = = t°. The contribution of sec- 
ond parameter may be observed on the right hand side of 
equation (23). For f P =Q and n = 1, equation (23) reduces 
(except for a minus sign) to the equation derived by WS (cf. 
eq.(23) in WS). When /, < 1 near the centre of the disk, 
we find 



(t)' 



l +Jt UJ ■ m 

Toward the centre, A increases which shows that gas surface 
density decreases (r varies inversely with A). Large values 
of fp for diffuse clouds again guarantees the depletion of 
gas in the centre. For outer parts of disk, f g ~ 1, and af- 
ter expanding various terms in equation (23) and neglecting 
higher order terms, we get the following 

V 



(1 + /p)/5 



+ 1 



(25) 

<§C 1, the gas density 



In view of large values of r and /, 
scales as the initial one. 

We now proceed to obtain critical column densities 
based on n and that on new parameter Qa- We write 

KV ef f 



O'er, A 



■KG 

2V2Av eff 
nG 



(26) 



Assume a rotation curve of the form V cx r M (n — for flat 
curve) . One gets 

Ocr. A 1 



/' 



(27) 



It is found that for fi = 0, both densities agree. But for large 
H (i.e. departures from flatness), o cr ,A becomes smaller than 
o C r,K- For example, for M 33, u = 0.3 (Newton 1980) which 
yields 



O'er, A = 0.61<7cr,K • 



(28) 



Observations of o g (Wilson, Scoville & Rice 1991) for this 
galaxy are better explained if one takes o cr ,A as threshold 
density rather than o cr , K (see also Elmegreen 1993). Thus, 
Qa emerges as a better parameter relating to disk instabil- 
ities than Toomre Q-parameter for star formation. This is 
also supported by the ratio of the two threshold densities 
(see e.g. Table 2). 

Table 2 shows that for highly non-linear region of rota- 
tion velocity, threshold density based on Qa is lowered (rela- 
tively) favouring instability for star formation. On the other 
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Table 2. Variation of ratio of threshold densities with index fi 







0.005 


0.99 


0.05 


0.93 


0.10 


0.86 


0.15 


0.80 


0.20 


0.73 


0.30 


0.61 


0.40 


0.51 


0.50 


0.41 


0.60 


0.32 


0.70 


0.23 


0.80 


0.12 


0.90 


0.07 



1.9 




< 

o' 



0.8 



003 








Table 3. Radial variation of Q 


a/Q for the Galaxy 




r 


A 


-B 




Qa/Q 


(kpc) 


(. 1Q -16.5 S -1) (. 1Q 


-16.5 S -1) 


(.10-16.5 S -1) 




1 


105 


62 


203.5 


1.5 


2 


.30 


55 


136.8 


0.6 


3 


20.9 


44.5 


107.9 


0.5 


4 


19.7 


34.1 


85.7 


0.6 


5 


19.1 


26.1 


68.7 


0.8 


6 


18.2 


20.1 


55.5 


0.9 


7 


17.2 


15.6 


45.2 


1.1 


10 


13.8 


7.9 


26.2 


1.5 


12 


11.5 


5.5 


19.3 


1.7 


11 


9.6 


4.3 


15.5 


1.7 


16 


7.9 


3.6 


12.9 


1.7 


18 


6.5 


3.3 


11.4 


1.6 


20 


5.44 


3.11 


10.3 


1.5 


30 


2.91 


2.40 


7.1 


1.2 


50 


1.59 


1.53 


4.4 


1.0 


75 


1.06 


1.00 


2.9 


1.0 



0.3 



0.0 



20.0 



40.0 
d/kpc 



60.0 



hand, threshold density based on Q-parameter is higher 
(by about one order) in this region. This shows that Q- 
parameter is relatively less efficient to favour star formation. 
We, therefore, conclude that in the non-linear regime of ro- 
tation velocity curve, the Q-parameter is approximately 10% 
less effective than Qa parameter for triggering the process. 

We have computed the ratio Qa/Q (see e.g. Table 3, 
data taken from Einasto 1979) at various radial distances of 
the Galaxy. B denotes the second Oort constant. A plot of 
Qa/Q with radial distances from the centre is shown in Fig. 
1. It is found that Qa depicts almost similar behaviour as the 
observed Q-distribution (cf. Fig. 6 in WS) remarkably for the 
range 1 kpc to 15 kpc and stays at Qa/Q > 1 beyond 30 kpc. 
Thus, Qa and Q both parameters agree beyond 30 kpc, i.e., 
in the flat rotation curve region. We obtain same result from 
the data shown in Table 2. In fact, for a disk radius below 
30 kpc deviations in the two parameters become significant 
which shows the relative merit of Qa parameter over Q- 
parameter to keep track of the physical process like star 
formation and other nuclear activity as well. 



Figure 1. Radial variation of Qa/Q for the Galaxy. We have 
taken data from Einasto (1979) and have expressed Qa parameter 
as in Elmegreen (1993). 



3.1 COMPARISON OF OBSERVATIONS FOR 
THE GALAXY 

We assume constant IMF in the solar neighbourhood (Miller 
& Scalo 1979; Scalo 1986) and take the following input 
parameters: initial surface density ct^.q ~ a St Q + ct Si q ~ 
50M©pc" 2 (Kuijken & Gilmore 1989; Bahcall, Flymn & 
Gould 1992), G a ~ 10M e pc~ 2 (McKee 1990), f g ~ 0.2, 
/ ~ 0.05 (Elmegreen 1993), t = age of the Galaxy = 15 Gyr, 
a = 0.1 (Myers et al. 1986), Oort shear constant A — 15 
km s' 1 kpc -1 (Kerr & Lynden-Bell 1986), 5 = 0.3 (Miller 
& Scalo 1979; Scalo 1986). For n = 1, we get timescale 
of star formation as r = 0.38 Gyr. We calcualte fd,o us- 
ing equation (23) as fd,@ ~ 0.10. After substituting these 
values into equation (11), we get star formation rate as 
R = 3.8MQpc~ 2 Gyr~ 1 . This is in agreement with Scalo 
(1986) who infers R ~ (1 - 4)Af0pc" 2 Gyr _1 within an un- 
certainly factor of about 3. For n = 2, we get r = 0.12 Gyr, 
fd,Q ~ 0.05 and star formation rate R = 6.0MQpc~ 2 Gyr~ 1 . 

We find that our model with n = 1 provides star for- 
mation rate which is in good agreement with inferred rate 
in the solar neighbourhood. It is to be remarked that our 
models are sensitive enough to efficiency a introduced in 
equation (1) which however is determined by star formation 
timescale r. We note that even if the efficiency drops by 
10% (i.e. when the value of a becomes of the order of 0.01) 
model with n = 2 gives same star formation rate as model 
with n = 1 and a = 0.1. Parametric freedom for a and fp 
even when Qa > 1 (i.e. non-gravitational instabilities are 
dominant) provide a general star formation scenario. Our 
model thus presents a generalisation of WS model with a 
dependence of star formation rate on Oort shear constant 
A. In contrast to WS, we find continuous (in the sense of 
Q-values) star formation rate obeying a similar but however 
different criterion (i.e. Qa < 1) of gravitational instability 
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for gaseous disks. In fact, competitive nature of two terms 
in equation (15) helps one to visualize the essence of conti- 
nuity in the star formation process. We discuss the scenario 
in more details in Sect. 3.3. 

3.2 TIMESCALE OF GAS DEPLETION 

For a particular n, we get from equations (1) and (16) as 

= -a(l - 5) n r- n dt n . (29) 
Integrate equation (29) to obtain 

In iTj = —a(l — S) n T n t n + constant. (30) 

At t = 0, <7g(r,i) = <Tg{r,0) which yields 

<(r, t) = <(r, 0) exp [-a(l - *)"T-*t n ] (31) 

(see also Lynden-Bell 1975; Glisten & Mezger 1983). Now 
we can write e-folding time as 

*3 = a (l_i)« T -™ ' (32) 

For our input parameters, the depletion time td for model 
n = 1 is td — 5.4 Gyr and for model n = 2, td — 0.54 Gyr. 
Model with n = 2 has 10% depletion time as compared to 
n — 1. For an age of 15 Gyr of the Galaxy model with n = 1 
implies that present gas fraction is ~ 10% of its initial value 
assuming that there is little variation over the last 5 Gyr 
(Dopita 1985, 1987). 

3.3 THE /p-PARAMETER AND STAR 
FORMATION 

The /p-parameter introduced in equation (15) requires fur- 
ther analysis as regards the process of star formation. It is 
dimensionless and measures the fraction of diffuse clouds 
to self-gravitating clouds. Low values of fp (fp ~ 0.01) 
means that clouds are dense and self-gravitating. In this 
case, physics of star formation is largely determined by 
the first term in equation (15). If, however, fp ~ 100 as 
for example in the inner Galaxy where pressure becomes 
high (Elmegreen & Elmegreen 1987, Polk et al. 1988, see 
also Vogel, Kulkarni & Scoville 1988 for M 51), as a result 
diffuse molecular clouds collide and cool leading to large 
mass cloud formation. Nevertheless, this does not mean that 
such regions result into large star- forming clouds. In fact, 
gravitational instabilities are more efficient (as compared to 
diffuse cloud collision) to produce large mass star-forming 
clouds. But, then, in this case local energy dissipation oc- 
curs through diffuse cloud collisions (Elmegreen 1989). A 
major difficulty for gravitational instability triggered star 
formation appears when Qa and fp both are large. In this 
case, only thermal instability is responsible for switching on 
star formation process. Murray & Lin (1989) have stressed 
the domonating role of thermal instability over gravitational 
instability for a proto-globular cluster where fragmentation 
(into protostars) is initiated by the former. Low fp values 
may also result when pressure becomes low (i.e. in the outer 
spiral arms of galaxies where gravity is not significant to 
form large molecular clouds) and star formation proceeds 
via shear instability. This instability does not depend upon 



Table 4. Radial variation of surface density for the Galaxy 



d 

(kpc) 


log ai(M Q pc 2 ) 


log o g {M e pc 2 ) 


0.001 


5.5 




0.01 






0.1 


4.4 




1.0 


3.0 


2.0 


2.0 


3.0 


0.5 


3 


2.8 


0.7 


4 


2.7 


1.0 


5 


2.5 


1.0 


6 


2.4 


1.0 


7 


2.2 


1.0 


10 


1.7 


0.9 


12 


1.4 


0.7 


14 


1.0 


0.6 


16 


0.7 


0.5 


18 


0.4 


0.3 


20 


0.1 


0.0 



Qa- Still, Qa has to be relatively small to guarantee un- 
stable radial motion which in turn facilitates dense cloud 
formation. Hence, the new parameter fp appearing in equa- 
tion (15) also explains star formation scenario plausibly in 
the spiral arms and outer parts of galaxies. 

We suggest that two competitive parameters {Qa and 
fp) may be understood to ultimately decide star formation 
(in view of observed continuity of the process from centre to 
arm for the Galaxy) as follows. We propose a relation 

fp = 10" 1 " qA for Q A < 1 (33) 
f P = 10 1+Qa for Qa > 1 • 

From equations (33), we find that Qa changes continuously 
from its value — (1 + log fp) to (— 1 + log/p) respectively 
from gravitational instability region (fp ~ 0.01) to non- 
gravitational instability region (fp ~ 100) and star forma- 
tion is supposed to occur from near the centre to outer spi- 
ral arms in contrast to earlier suggestion (Quirk 1972) who 
favours (necessarily) gravitational instability to trigger the 
process. 



4 VARIATION IN THE STAR FORMATION 
RATE 

We assume cloud mass density in the solar neighbourhood, 
a c ~ IOOMqjjc -2 , as constant (Larson 1981). We further 
assume / c ~ 0.01 at d = 1 kpc (since a g at 1 kpc is 
~ 100 Mq pc~ which makes f c = 1 yielding an infinite R) 
to keep star formation reasonably large in the model. Com- 
ponents of surface densities (data taken from Einasto 1979; 
Lacey & Fall 1985 and WS) are given in Table 4 and plotted 
as Figure 2 at various distances from the Galactic centre. 

We infer from Figure 2 that the Einasto model shows 
ct, ~ r -°- 8 dependence for r < 6 kpc and deviates for 
r > 6 kpc (see also e.g. Kundt 1990 for a variant in the 
Galactic mass distribution). Now, we aim to discuss the vari- 
ation of star formation rate normalised to that in the solar 
neighbourhood and therefore calculate R/Rq (see e.g. Table 
5) using data from Einasto (1979) and plot as Figure 3 at 
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Figure 2. Radial variation of total surface density ai(MQpc~ 2 ) 
and gas surface density tj g (MQpc~ 2 ) for the Galaxy. The solid 
line represents the data taken from Lacey & Fall (1985) and WS; 
the dotted line represents the data from Einasto (1979). 



Figure 3. Star formation rates normalized to its value in the solar 
neighbourhood. The data are based on Einasto (1979); Lacey & 
Fall (1985) and WS. 



Table 5. Radial variation of star formation rate for the Galaxy 



d 


A 


9% 


fc 


R/Rq 


(kpc) 


( . 10 -16.5 S -1) 


(M e pc- 2 ) 






1 


105 


1016.3 


0.01 


78.6 


2 


30 


851.1 


0.03 


56.5 


3 


20.9 


633.9 


0.05 


48.8 


4 


19.7 


452.9 


0.10 


66.0 


5 


19.1 


318.4 


0.11 


49.5 


6 


18.2 


222.3 


0.10 


29.9 


7 


17.2 


154.9 


0.10 


19.7 


10 


13.8 


50.9 


0.07 


3.6 


12 


11.5 


23.6 


0.05 


1.0 


14 


9.6 


10.7 


0.04 


0.3 


16 


7.9 


4.9 


0.03 


0.1 


18 


6.5 


2.3 


0.02 


0.02 


20 


5.44 


1.2 


0.01 


0.005 



various Galactocentric distances. From Figures 2 and 3 we 
infer that star formation rate varies like gas component of 
surface density. 

A minimum in a g occurs at ~ 3 kpc where we also 
observe a minimum in the star formation rate. Thereafter, 
a g increases again and reaches a maximum at ~ 4 kpc where 
we observe corresponding increase and maximum in R/Rq. 
Our model may be applied to a general dependence of a g on 
the exponent n. A given value of n determines the timescale 
of star formation (r) which however yields the corresponding 
R. For n — 1, our model agrees with WS model but we 
obtain larger Schmidt exponent (see e.g. equation (12)). 

The star formation rates inferred from (i) pulsar data 
(Lyne, Manchester & Taylor 1985), (ii) from observations of 



supernova remnants (Guibert, Lequeux & Viallefond 1978), 
and (iii) Lyman-continuum photon observations from H II 
regions (Glisten & Mezger 1983) are consistent from our 
model at all radial distances. For example, higher rate of star 
formation traced by Lyman-continuum near 4 kpc agrees 
with our model calculations. This is evident by the maxi- 
mum in Figure 3 at 4 kpc from the Galactic centre. In view 
of comments (Wyse & Silk 1989) regarding higher star for- 
mation rates of Glisten & Mezger than those given by Scalo 
(1988) (viz. these estimates may be higher by an order of 
magnitude) and also the fact that it does not match with 
star formation profile obtained by other techniques (Rana 
& Wilkinson 1986), our values are apparently better tuned. 

We assume r = 0.45 Gyr. Using parameters as decribed 
in Section 3.1, we calculate efficiency a of star formation as a 
function of distance from the galactic centre. It is interesting 
to observe that a changes in the solar neighbourhood. Small 
values of a at 1 kpc may be understood to arise because 
of shear instability which removes growth of perturbations. 
Star formation can proceed if the self-gravitational collapse 
time becomes shorter than the shear time (~ 0.01 Gyr). 
However, relatively large value of a out to 10 kpc does not 
lead to large star formation rate R/Rq (see Figure 3) due to 
paucity of gas. In fact, density <j, drops below the observed 
value (Wilson, Scoville & Rice 1991) of critical density at 
14 kpc where we expect turn-off of star formation due to 
gravitational instability. It is also supported by significant 
depletion of gas at this distance (see e.g. Figure 2). The 
striking feature of our result is that a changes in the solar 
neighbourhood (an efficiency gradient ~ 0.0057fcpc~ 1 ) by an 
amount ~ 0.02. It seems thus natural to think that the effi- 
ciency gradient is responsible for radial abundance gradients 
which are reported in many disk galaxies (Edmunds & Pagel 
1984; Diaz & Tosi 1984; Tosi & Diaz 1985). The fact that 
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metallicity gradients may be due to changes in efficiency of 
star formation was suggested previously by Lacey & Fall 
(1985). We aim to confirm this suggestion from our calcula- 
tions too. There is hardly a need to invoke radial flows (see 
discussion in Scalo 1988) in this scenario. 

4.1 METALLICITY GRADIENT VS 
EFFICIENCY GRADIENT 

Following Pagel & Patchett (1975) (see also e.g. Pagel & 
Edmunds 1981) model of chemical evolution of the Galaxy 
in the solar neighbourhood, we define £, a mass ratio in the 
form of long-lived stars, and p as the yield of heavy elements 
which represents mass ejected per unit mass of long-lived 
stars (cf. Searle & Sargent 1972; Talbot & Arnett 1973a). 
For our model o s = £i7i, a g = (1 — £)ct;, a s = £/(l — £)a g 
and 



do s 
lit 



-£ dt 



(34) 



Equation (34) gives 

d n , "(1-0 /o K x 

-(lna 9 ) = — . (35) 

Integration of equation (35) yields 

In o„ = — _|_ cons i an i^ <2,q\ 

At t = 0, a g = Oi hence constant = In/Ji. Thus, equation 
(36) takes the form 



"(i-O t 



^ = (l-0=exp 

Oi 

Therefore, 



(37) 



(38) 



Metallicity Z is expressed as (Pagel & Patchett 1975) 



Z =pln 



P«(l ~ « 



(39) 



We see that r is now a function of time and is given by 
equation (38). Time evolution of Z may be written as 



(40) 



dZ _P _ P§_ 

~dt ~ ~ ~ a(l-0[ln(l-0] _1 * ' 
From equation (39) we infer that Z varies linearly both with 
time and efficiency of star formation. We assume that £ ~ 
0.8 (Talbot & Arnett 1973a) and p ~ 0.7Z Q (Wang & Silk 
1993) to calculate Z in the solar neighbourhood. For an 
efficiency of a ~ 0.07 we find Z ~ 1.13Zq for solar age. 
This is in agreement with the plot of metallicity in the solar 
neighbourhood by Wyse & Silk (1989, cf. Figure 2b). The 
present model thus provides time evolution of metallicity 
which however depends upon efficiency of star formation. 

For an efficiency run of a ~ 0.07,0.08,0.09,0.10, we 
find Z/Zq ~ 1.23,1.17,1.13,1.09 respectively which is in- 
depedent of Galactic age provided of course the parameters 
p and £ do not change with time. In other words, disk ag- 
ing alters r(t) such that t/r(t) remains constant (for a fixed 
a) and hence there occurs no change in Z/Zq. Our calcu- 
lations show that metallicity decreases with increase of a 



at a given age. In the solar neighbourhood, this may be 
understood due to paucity of gas favouring relatively low 
star formation at large distances (see Figure 3) and there- 
fore low metal production (see also Friel & Janesl993). We 
note that the observed run of metallicity of G-K dwarfs in 
our Galaxy is very sensitive to chemical composition of stars 
of same age (Tinsley 1975). Janes & McClure (1972) have 
suggested enhancement in the dispersion due to chemical in- 
homogeneities in the Galaxy (Talbot & Arnett 1973b). The 
structure of Galactic disk and the presence of population 
gradients are given in Ferrini et al. (1994). For a radial dis- 
tribution of abundances in galaxies one may refer to Molla et 
al. (1996). They have also discussed the chemical evolution 
of solar neighbourhood, see e.g. Pardi et al. (1995). 

However, as remarked by Tinsley (1975) that the ob- 
served dispersion (see also Hearnshaw 1972 for dispersion) in 
metallicity for stars of same age may result either partly due 
to chemical inhomogeneities (of interstellar medium) or due 
to causes altogether different, essentially favours this analy- 
sis. We find that metallicity dispersion for stars of same age 
may be due to variation of efficiency a through which differ- 
ent sample stars were processed. This confirms the assump- 
tion of Rana & Wilkinson (1986) that metallicity dispersion 
is due to stellar processing only. It is found that a depends 
upon star formation rate and also the gas component of sur- 
face density o g . We conclude that a predominantly deter- 
mines the observed dispersion and plays a key role towards 
metal enrichment or otherwise of the interstallar medium. 

At various disk ages (at a given radial distance), there 
occurs change in a which causes metallicity dispersion. We 
notice that a also changes at various distances from the 
Galactic centre which results in spatial metallicity gradi- 
ents. One immediately finds that apparent metallicity dis- 
persions with either age or distance depend upon a. The 
[O/H] vs age plot (Wyse & Silk 19879, cf. Figure 2d; see 
also Carlberg et al. 1985) shows hardly significant metal- 
licity gradient at all disk ages (cf. Friel & Janes 1993). We 
suggest that all sample stars might have followed evolution 
with almost similar efficiency. Thus, the important result of 
this analysis is the confirmation of the suggestion by Lacey 
& Fall (1985) and Richtler (1995) regarding metallicity gra- 
dients. For a comprehensive treatment of radial abundance 
gradients in spiral disks and age-metallicity relation in dif- 
ferent stellar populations one may refer to Edvardsson et 
al. (1993); Pagel (1994). An interesting modern analysis of 
kinematics and abundance distribution for our Galaxy has 
been given by Gilmore, Wyse & Kuijken (1989). Matteucci 
(1996) has reviewed exhaustively the evolution of the abun- 
dances of heavy elements in gas and stars (indicating ob- 
servational and theoretical constraints) in galaxies of differ- 
ent morphological types. After a similar work by Tinsley 
(1980), the article provide a good document of the progress 
in the understanding of the physical processes regulating the 
chemical evolution of galaxies. Formation and evolution of 
our Galaxy is also discussed. For a review on abundance ra- 
tios and Galactic chemical evolution see Mc William (1997). 
Chemical evolution of solar neighbourhood according to the 
standard infall model using SN II data are summarized in 
Thomas et al. (1998). 
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5 DISCUSSIONS 

We have studied that the suggestion by Elmegreen (1993) 
regarding star formation appears more robust than Q- 
criterion. It is because unless Q < 1, gravitational instability 
does not permit star formation. However, when Q > 1, the 
system becomes gravitationally stable and consequently star 
formation via large cloud formation is not feasible. A natural 
question which one would ask is: how does star formation 
proceed when Q enters the stable regime? This infact led 
to an alternative criterion for cloud formation (discussed in 
the text) leading to star formation as originally suggested by 
Elmegreen (1993). Accordingly, when magnetic field is taken 
into consideration, velocity dispersion changes and thus Q is 
pushed to the stable regime. At this stage, non-gravitational 
instabilities (e.g. thermal instability, shear instability) dom- 
inate over gravitational instability. We infer from Figure 1 
that dependence of Qa/Q on distance from the Galactic 
centre describes the relative merit of Q^-parameter over Q- 
parameter beyond 6 kpc. Observations of a g for M 33 are in 
better agreement with theory when one regards Qa as the 
gravitational instability parameter (see e.g. Wilson, Scoville 

6 Rice 1991). It is found that both Qa and Q parameters 
agree beyond 30 kpc. 

We have obtained generalised version of WS equation 
(see e.g. Sect. 2 equation (12) and Sect. 3 equation (23)) 
in the sense that (i) there is an additional non-zero term in 
equation (12) and (ii) in view of equation (15), one arrives at 
a natural rescue from cut-off criterion for star formation. We 
also show that the nominal Schmidt exponent n s is given by 
2 < n s < 3 in our model. We suggest a general equation (e.g. 
equation (15)) for the star formation rate consisting of two 
terms: the first term dominates when Qa < 1 and fp <JC 1; 
the second term dominates when Qa > 1 and fp <C 1. 
Apparently, the relative competence of either of these terms 
determines star formation scenario as discussed in Sect. 3 
at all radial distances. Virtually, Q < 1 (or Qa < 1) is 
not an absolute criterion for star formation. For our model 
with n = 1, we get star formation rates which are in good 
agreement with values inferred by Scalo (1986). We find that 
our models are sensitive enough to efficiency a and timescale 
t of star formation. A given exponent n determines r which 
however yields the corresponding star formation rate R. 

We suggest that essentially the efficiency gradient is the 
cause for radial abundance gradients which are reported in 
many disk galaxies. Under the approximation of closed box 
model, we have derived time evolution of r(t) and also the 
metallicity Z(t). Both r(t) and Z(t) are functions of a, p 
(the yield of heavy element) and mass ratio £. We notice 
hardly any metallicity change as disk ages which however 
reflects that stellar processing occurs at a fixed a. Metallic- 
ity dispersion for stars of same age may be caused due to 
variation in a. We conclude that a is predominantly respon- 
sible for metallicity dispersion and also for metal enrichment 
of interstellar medium. A simple model as discussed above 
provides some important characteristics of our Galactic disk. 
Although, as suggested by Tinsley (1980), the star forma- 
tion is a complicated function of several physical parameters, 
e.g. gas density, gas sound speed, shock frequency, shock 
strength, gas rotation, shear constant A, magnetic field, gas 
metal abundance, and background star density. It is how- 
ever difficult to predict the actual dependence of R on these 



parameters. One therefore studies some form of R and its 
consequent effect on chemical and photometric evolution. 
Finally, the model predictions are compared with observa- 
tions. 

We note that star formation rate was probably higher 
in the central part of the disk of our Galaxy at an early 
epoch of evolution. It is to be remarked that hydrodynam- 
ical simulations of the formation and evolution of a galaxy 
may be performed incorporating our model formulation of 
star formation rate and metallicity. Model predictions when 
compared with observations of other galaxies would speak 
of its robustness and proof. 
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